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Abstract. We report on observations of the recently detected (Kuzmin et al. I2004t giant pulses (CPs) from the 
pulsar PSR B0031-07 at 40 and 111 MHz. At 40 MHz the peak flux density of the strongest pulse is 1 100 Jy, 
which is 400 times as high as the peak flux density of the average pulse (AP). A pulse whose observed peak flux 
exceeded the peak of the AP by more than a factor of 200 is encountered approximately once in 800 observed 
periods. Peak flux density of the CPs compared to the AP peak flux density <S'^oak/'5']^e^k has roughly a power-law 
distribution with a slope of —4.5. CPs at 40 MHz are essentially stronger than those ones at 111 MHz. This excess 
is approximately in inverse proportion to the frequency ratio. The giant pulses are much narrower than the AP, 
and cluster in two narrow regions of the AP near the peaks of the two components of the AP. Some of the CPs 
emit at both phases and are double. The separation of the double CP emission regions depends on frequency. 
Similarly to the frequency dependence of the width of the AP, it is less at 111 MHz than at 40 MHz. This suggests 
that CPs are emitted from the same region of the magnetosphere as the AP, that is in a hollow cone over the polar 
cap instead of the light cylinder region. PSR B0031-07 as well as the previously detected PSR B1112+50 are the 
first pulsars with CPs that do not have a high magnetic field at the light cylinder. One may suggest that there 
are two classes of CPs, one associated with high-energy emission from outer gaps, the other associated with polar 
radio emission. The CPs of PSR B0031-07 and PSR B1112-I-50 are of the second class. The dispersion measure 
DM is found to be 10.900 ± 0.003 pc cm"^. 
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1. Introduction 

Giant pulses (GPs) are short-duration burst-like sudden 
increases of intensity of individual radio pulses from pul- 
sars. This rare phenomenon has been detected previously 
only in five pulsars out of more than 1 500 known ones:: 
the Crab pulsar PSR B0531-h21 (StaeHn & Sutton 1^717), 
PSR B1112+50 (Ershov & Kuzmin 2003J, the miUisecond 
pulsars PSR B1937-f-21 (Wolsz czan e t ai. IlM^ and PSR 
B1821~24 (Romani & Johnston l2001|l . and the extragalac- 
tic pulsar PSR B0540-69 (Johnston & Romani 2003). 

The intensities of the GPs are extremely high. The 
brightness temperature of GPs of the millisecond pulsar 
PSR B1937-h21 is as high as Tb w lO^^ K (Kuzmin & 
LosovskyEOni, and Tb > lO^^ K for GPs of the Crab 
pulsar (Kostyuk et al. 12003 Hankins et al. l2003p . 

Four of these pulsars belong to the group of pulsars 
with the strongest magnetic field at the light cylinder. 
Pulsar PSR B1112+50 is the first exception from this 
group, being a pulsar that exhibits the characteristic fea- 
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tures of GPs but has an ordinary value of the magnetic 
field on the light cylinder. Kuzmin et al. H2004|l recently re- 
ported the detection of giant pulses from another ordinary 
magnetic field pulsar, PSR B0031-07. They find that at 
111 MHz the peak flux density of the strongest CP is 530 
Jy, which is 120 times as high as than the peak flux density 
of the average profile (AP). A pulse whose observed peak 
fluxes exceed the peak of the AP by more than a factor 
of 50 is encountered approximately once in 250 observed 
periods. Giant pulses cluster in a narrow region near the 
peak of the first component of the AP. 

In this paper we report observations of GPs of pul- 
sar PSR B0031-07 at the second frequency of 40 MHz, 
additional observations at 111 MHz and a comparative 
analysis of the two-frequency data. 

2. Observations and data reduction 

Observations were carried out at two frequencies, 40 and 
111 MHz. 

At 40 MHz observations were performed from August 
28 through October 13, 2003 with the DKR Radio 
Telescope at Pushchino Radio Astronomy Observatory of 
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the Lebedev Physical Institute. The DKR telescope has 
an effective area of about 8000 square meters. One lin- 
ear polarization was received. We used a 128-channel re- 
ceiver with channel bandwidth 1.25 kHz. The frequency of 
the first (highest frequency) channel was 40.982 MHz, the 
sampling interval was 0.819 ms, and the time constant was 
T = 1 ms. The duration of each observation was about 16 
min (994 pulsar periods). A total of 16 observations con- 
taining 15 904 pulsar periods was carried out. 

From June 21 through September 22, 2003 we per- 
formed an additional 16 observations at 111 MHz and 
12 simultaneous observations at 111 and 40 MHz. At 
111 MHz we used the Large Phased Array (BSA) Radio 
Telescope at Pushchino Radio Astronomy Observatory 
with an effective area of about 20 000 square meters. One 
linear polarization was received. At both frequencies we 
used two 64-channel receivers with channel bandwidth 20 
kHz. The sampling interval was 2.560 ms and the time 
constant was r = 3 ms. The durations of the observations 
were about 3 min (205 pulsar periods) and 9 min (600 
pulsar periods) at 111 and 40 MHz respectively. 

All observations were time-referenced to the 
Observatory rubidium master clock, which in turn 
was monitored against the National Time standard via 
TV timing signals. Pulsar ephemerides were obtained 
from the catalogue by Taylor et al. H1995|l . 

During the off-line data reduction the signal records 
were cleaned of radio interferences. Subsequently the 
inter-channel dispersion delays imposed by an interstellar 
medium were removed. Each observing period was ana- 
lyzed for pulses with amplitude exceeding a preset level 
and its amplitude, pulse width and phase were derived. 



3. Results 

At 40 MHz the 212 pulses whose observed peak fluxes den- 
sity exceeded the peak flux density of the average profile 
(AP) (which is equivalent to the average pulse) by more 
than a factor of 100 were selected and analyzed. 19 of 
them exceeded the peak of the AP by more than a factor 
of 200. A pulse whose observed peak fluxes exceeded the 
peak of the AP by more than a factor of 100 is encountered 
approximately once in 75 observed periods. 

Figure^ shows the strongest observed GP (bold line) 
together with the AP for all 16 days of observations. The 
observed peak flux density of this GP exceeds the peak 
flux density of the AP by factor of 400. The plot of the AP 
is presented on a 200 times larger scale and flux densities 
of the observed GP and AP are shown separately on the 
left and right sides of the "y" scales. 

The value of the GP flux density was determined rel- 
ative to the known AP flux density as 




cGP 
peak 



cAP 
peak 



X (/GP//AP) 



where /GP//AP ^j^g ^.j^^Jq ^j^^ Qp ^-^ pgj^j^ 
intensity and S'^^^k the peak flux density of AP 
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Fig. 1. a) The strongest observed GP at 40 MHz (bold 
line) together with the AP (thin hue). The plot of the AP 
is presented on a 200 times larger scale and flux densities 
of the observed GP and AP are shown separately on the 
left and right sides of the "y" scales, b) The phases of the 
observed GPs relative to the center of the AP determined 
as the mid-point between 0.1 intensity levels. 



cAP _ cAP /, 
'^pcak ~ '^nioan/ "'form j 

where S^oan is the flux density averaged over a pulsar 
period, feform is the duty cycle of the pulsar. 

The value of the AP flux density 5",^,^^,, = 490 mJy 
was obtained from Izvekova et al. (1981). The reduction 
factor is feform = 0.18, so the AP peak flux density is 2.7 Jy. 
The strongest observed GP has /GP//AP ^ 490. Then, the 
observed GP peak flux density is S^^^^ = 2.7 x 400 = 1 100 
Jy. The intrinsic flne structure of the pulses is masked 
by interstellar scattering. In our measurements of the 19 
strongest GPs at 40.9 MHz we derived r^c = 3.1 ± 0.4 
ms. Scattering increases the pulse width and decreases the 
pulse amplitude. Thus the intrinsic GP peak flux density 
is 5pGP^ > 1 100 Jy 

In Fig. |2t we show the cumulative distribution of the 
ratio of the observed GP peak flux density to the AP peak 
flux density ■S'^oak/'^'^cak '^^ pulses that we have selected 
and analyzed. In a detectable interval of ■S'^oak/'^'^eak from 
100 to 400, the histogram has roughly a power-law distri- 
bution with a slope of -4.5. Since for S^^ak/'^'peak < 100 
the pulse signal-to- noise ratio is less than 5, at lower in- 
tensities the observed distribution is masked by noise. 

The distribution for 111 MHz, which we have obtained 
from the Kuzmin et al. 1)2004(1 data and additional obser- 
vations, is shown in Fig. Eb. In a detectable interval of 
'S'^cak/'^pcak from 30 to 100, the histogram has roughly 
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Fig. 2. a) The cumulative distribution of the observed GP 
peak flux density as related to the AP peak flux density 
for 40 MHz, b) The same for f If MHz. 



a power-law distribution with a slope of —4.8. Since for 
'^^ak/'S'peak < 20 the pulsc signal-to-noisc ratio is less 
than 5, at lower intensities the observed distribution is 
masked by noise. 

The distributions in Fig. |21 demonstrate, that GPs at 
40 MHz are essentially stronger than those at 111 MHz. 
At 40 MHz one in every thousand pulses exceeds the AP 
by more than a factor of 170, but at 111 MHz only by 
60. The strongest observed GP at 40 MHz exceeds the 
peak flux density of the AP by a factor of 400, but at 111 
MHz only by a factor of 120. In a quantitative sense this 
gain is approximately in inverse proportion to the ratio of 
frequencies. 

From Fig. |2] one can see that the rate of GPs depends 
on the frequency. At 40 MHz a pulse whose observed peak 
flux exceeds the peak of AP by more than a factor 100 is 
encountered once in 85 observed periods, but at 111 MHz 
only once in 11 000. This may be a reason why there are 
no claims of GPs from this pulsar at higher frequencies. 

Our distributions cover a smaller number of observed 
pulsar periods (by two orders of magnitude) compared to 
such histograms for previously known short period pul- 
sars. Therefore, we consider the slope index as a tentative 
one which needs to be refined. 



The mean observed width^ of 19 GPs with 
^pTak/^p^cak > 200 is u^st^ = 6 ± 2 ms. The width of 
the AP^ for all 12 days of observations is = 180 ms. 
Thus, at 40 MHz GPs are narrower than the AP by about 
a factor of 30. 

The positions of the GPs are stable within the AP and 
they cluster in a narrow window near the maximum of the 
second component of the AP. Figure ^Jj shows the phases 
$ of the GPs which exceeded the AP by factor of 100. 
The average phase difference between the GPs and the 
center of the AP is $40 = 40 ms. GPs cluster in a narrow 
phase window A<i> = ±10 ms. The clustering is tighter for 
stronger GPs. 

The GPs are not identified with subpulses, since con- 
trary to the stable phase position of GPs, the phase po- 
sition of subpulses drift across the AP. GPs are also nar- 
rower than the subpulses of this pulsar w^q = 15 ms at 60 
and 102 MHz (Izvekova et al. lT^^ . 

The GP brightness temperature is 

Tb = 5AV2fcrj . 

Here A is the radio wavelength, k is the Boltzmann 
constant, and Q is the solid angle of the radio emission 
region. Adopting ~ (l/d)'^, where I is the size of the 
radio emission region and d the distance to the pulsar, 
and further adopting I < c x W50, where c is the speed of 
light and the distance to the pulsar d — 0.68 kpc (Taylor 
et al. 19a5i, we obtain Tb > 10^^ K. 

The most unexpected and interesting feature of our 
observations at 40 MHz is the non-coincidence of the phase 
positions of the GPs at 40 and 111 MHz. In Fig. El we 
present the time-alignment of GPs and APs at 40 and 111 
MHz (Kuzmin et al. lW)^ . One can see that the GP at 40 
MHz is separated from the GP at 111 MHz by about 60 
ms and that they are located in two different parts of the 
AP. At 111 MHz the GP is located near the maximum of 
the first component of the AP, whereas at 40 MHz the GP 
is located near the maximum of the second component. 

To verify this peculiarity, we have performed a con- 
trol determination of the dispersion measure DM. For 
this we varied the value of DM by a process of re- 
moval of the inter-channel dispersion delays and searched 
their optimal value, which provides the maximum am- 
plitude of individual pulses. We processed 440 individ- 
ual pulses at 111 MHz and 337 at 40 MHz. The mean 
value of dispersion measure was obtained to be DM ~ 
10.900 ± 0.003 pc cm^"^. This value is consistent with 
the catalog value 10.89 ± 0.01 pc cm~^ (Taylor et al. 
I1995|) and confirms the validity of the non-coincidence of 
the phase positions of the GPs at 40 and 111 MHz. 

One may suggest two explanations of such behavior. 
The first is that the emission region of GPs moves with 



^ The width is not corrected for scatter and dispersion broad- 
ening 

^ The width of the AP was measured at half the maximum 
height of the leading and trailing components. 
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Fig. 3. a) The strongest observed GP at 40 MHz (bold 
line), together with the AP (thin line), b) The strongest 
observed GP at 111 MHz (bold line), together with the 
AP (thin line). Plots of GPs and APs are presented on 
different scales, shown on the left and right sides of the 
"y" axis. 



frequency. The second proposes that there are two sepa- 
rate phase emission regions for the GPs. In this case, one 
may expect that some of the GPs will emit at both phases 
and are double. 

To clarify this point, we undertook a more detailed 
inspection of the GPs. Our first-line reduction program 
searched only for a single GP per pulsar period, whose 
amplitude exceeded a preset level. However, one may ex- 
pect that, in the case of two emission regions, some GPs 
have two components, where the first component is ob- 
served mainly at 111 MHz and the second at 40 MHz. 
Therefore we performed an addition search of the observed 
GPs for two-component structure. For this search we used 
a 20 kHz channel bandwidth at both frequencies. The an- 
swer is yes: — some GPs, both at 40 and 111 MHz, have 
a two-component structure. We detected 12 double GPs 
whose amplitude exceeded the AP by a factor of 50 at 111 



Phase, ms 

Fig. 4. a) The double GP (bold hne) observed at 111 
MHz together with the AP (thin line), b) The double GP 
(bold line), observed at 40 MHz. together with the AP 
(thin line). 



MHz and 5 double GPs with amplitude exceeding AP by 
a factor of 20 at 40 MHz. 

An example of two-component GPs is shown in Fig. 01 
Since at 40 MHz double pulses were searched in a 20 kHz 
channel bandwidth, the dispersion pulse broadening of 30 
ms reduces the strength of the GP and distorts the ob- 
served pulse shape compared to other GPs, which were 
observed in 1.25 kHz channel bandwidth. 

The phases of these double GPs, as well as single GPs 
from our observations and the Kuzmin et al. (^200111 data 
are plotted in Fig. [S] together with APs. 

One can see that the phases of the trailing component 
of double GPs at 40 MHz (Fig. Et) fit the phases of the 
single GPs at this frequency (Fig. [Sja) . In a similar man- 
ner, the phases of the leading component of double GPs 
at 111 MHz (Fig. |5ji) fit the phases of the single GPs at 
111 MHz (Fig.lSjj). This supports the concept of two emis- 
sion regions of GPs. The first dominates at 111 MHz, the 
second dominates at 40 MHz and both form double GPs. 
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Fig. 5. a) Integrated profile at 40 MHz, b) Phases of the 
observed single GPs at 40 MHz, c) Phases of the observed 
double GPs at 40 MHz, d) Phases of the observed double 
GPs at 111 MHz, e) Phases of the observed single GPs at 
111 MHz, f) Integrated profile at 111 MHz. 



Twelve observation sessions, containing 2 460 pulsar 
periods, were performed simultaneously at 40 and 111 
MHz. We also used the 20 kHz channel bandwidth at both 
frequencies. Forty pulses at 40 MHz and 45 pulses at 111 
MHz with peak intensities exceeding the peak intensity of 
the AP by more than a factor of 20 were selected and ana- 
lyzed. Only two among these pulses were found in the same 
pulse phase inside the same pulsar period and identified 
as simultaneous. The probability that there are acciden- 
tal coincidences of GPs simultaneously at two frequencies 
is equal to P40 * Pin = 40/2460 * 45/2460 « 10"^; this 
is significantly lower than the observed detection rate of 
simultaneous GPs of 2/2460 w lO'^. 

The measured spectral indices between frequencies 
of 40 and 111 MHz for these 2 simultaneous GPs are 
= —0.6 and —1.1, corresponding to the low-frequency 
turn-over region of this pulsar (the spectral index of the 



AP is a = —0.2). However, only a small number of si- 
multaneous GPs was observed. This possibility therefore 
needs further study. 



4. Discussion 

The GPs that we detected from PSR B0031-07 exhibit 
all characteristic features of the classical GPs from PSR 
B0531-h21 and PSR B1937-I-21. 

The peak intensities of the GPs exceed the peak in- 
tensity of the AP by more than a factor of 400. The his- 
tograms of the flux density have a power- law distribution. 
The GPs are much narrower than the AP and their phases 
are stable inside the integrated profile. 

Alongside these similarities, the GPs of the pulsar PSR 
B0031-07 have noticeable differences. Previously known 
pulsars with GPs belong to the group of pulsars with 
very high magnetic field on the light cylinder. Pulsar 
PSR B0031-07, as wefl as the previously detected PSR 
B1112+50 (Ershov & Kuzmin are the first ones 

with ordinary magnetic fields on the light cylinder. 

The most unexpected and interesting result is the clus- 
tering of the GPs in two different regions. This indicates 
that there are two emission regions of GPs. The separa- 
tion of these regions at 40 MHz is larger than at 111 MHz; 
this is corresponds to the increase in the width of the AP, 
which is interpreted as a divergence of the magnetic field 
lines in the hollow cone model of pulsar radio emission. 
This suggests that the GPs from this pulsar, and possibly 
from PSR Bl 112-1-50, originate in the same region as the 
AP, that is in a hollow cone over the polar cap instead of 
in the light cylinder region. 

One may suggest that there are two classes of GPs, one 
associated with high-energy emission from the outer gaps, 
the other associated with polar radio emission. The GPs 
of PSR B0031-07 and PSR B1112+50 are of the second 
class. 

5. Conclusions 

The detection of giant pulses from pulsar PSR B0031-07 
has been confirmed. The peak flux density of the strongest 
pulse is 1 100 Jy, which is 400 times as high as than the 
peak flux density of the average profile. 

The GPs are much narrower than the average profile 
and cluster in two narrow regions of the average profile. 
The separation of the two GP emission regions decreases 
with frequency. This suggests that the GPs from this pul- 
sar are emitted from the same region as the AP, that is 
in the hollow polar cone instead of in the light cylinder 
region. 

Pulsar PSR B0031-07, as well as the previously de- 
tected PSR B1112-I-50, are the first pulsars with GPs that 
do not have a high magnetic field at the light cylinder, such 
as previously known pulsars with GPs have. They may be 
a separate class of GPs. 

The dispersion measure DM is found to be 10.900 ± 
0.003 pc cm-3. 
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